TRACE ASYMPTOTICS FOR FRACTIONAL SCHRODINGER OPERATORS 



LUIS ACUNA VALVERDE 



Abstract. This paper proves an analogue of a result of Bafiuelos and Sa Barreto [6] on the 
asymptotic expansion for the trace of Schrodinger operators on R d when the Laplacian A, which 
is the generator of the Brownian motion, is replaced by the non-local integral operator A™/ 2 , 
< a < 2, which is the generator of the symmetric stable process of order a. These results also 
extend recent results of Bafiuelos and Yildirim [.'•>] where the first two coefficients for A a / 2 are 
computed. Some extensions to Schrodinger operators arising from relativistic stable and mixed 
stable processes are obtained. 



Heat asymptotic results have been widely used in areas of spectral theory and its applications to 
scattering theory, statistical and quantum mechanics and in several areas in geometry. We refer the 
reader to van den Berg [25] for the computation of the first two terms in the asymptotic expansion 
of the trace of the heat kernel of the Schrodinger operator — A + V under Holder continuity of 
the potential and to Bahuelos and Sa Barreto [G] for a more general computation with an explicit 
formula for all the coefficients for potentials V € iS(R d ), the class of rapidly decaying functions 
at infinity, and for applications to scattering theory. For applications in statistical mechanics 
and quantum theory, we refer the reader to the articles of Lieb [18] and Penrose and Stell [20] 
about the second viral coefficient of a hard-sphere gas at low temperature and sticky spheres, 
respectively. Heat trace asymptotic for the Laplacian have been of interest for many years for 



Supported in part by NSF Grant # 0603701-DMS under PI Rodrigo Bafiuelos. 



Contents 



1. Introduction. 

2. Stable Processes and subordinator. 

3. Heat trace in terms of Fourier transform. 

4. Boundedness of the (J+l)-th term. 

5. Heat trace computation by means of subordinators. 

6. Bounds for remainders and coefficients. 

7. An improvement for dimension d = 1,2,3. 

8. Proof of Theorem 1.1. 

9. Explicit form of some coefficients. 

10. Proof of Corollary 1.1. 

11. Explicit expansion fora = 2/fc, fc>2 integer and a close to 2 

12. Extension to a-relativistic processes. 

13. Extension to mixed-stable processes. 
References 



1 
7 
8 
11 
12 
15 
17 
19 
23 
2G 
27 
31 
33 
36 



1. Introduction. 



2 



LUIS ACUNA VALVERDE 



domains in Euclidean space M. d and on manifolds where the coefficients reveal many geometric 
quantities such as volume, surface area, convexity, number of holes, etc. For more on this large 
literature as well as some historical perspective, we refer the reader to Arndt, Nittka, Peter and 
Steiner, [ , pp 1-71], Bahuelos, Kulczycki and Siudeja [4, 5], Datchev and Hezari's [12], Donelly 
[14], McKean and Moerbeke [19], and Colin De Verdiere [11]. 

Let H2 — — A and Hy = — A + V, V € 5(K d ). In [6], the existence of an asymptotic expansion 
of the trace of the operator e~ tHv — e~ tH2 , as t 4- 0, is proved. To make the connection to the 
fractional Laplacian more clear, let us denote the heat kernel for —A by 



(47rf) 



d/2 



(47rf) 



d/2 ' 



so that 

Set 
(1.1) 

Throughout the paper we use the notation f(t) = 0(g(t)), as t I 0, to mean that there exist 
constants C and 5 such that \f{t)\ < C\g(t)\, for < t < S. 

With this notation the result in [6] can be stated as follows. For any integer J > 1, 

,7 



(1.2) 

as t 4- 0, with 

ci(V) 



Ij = {A = (Ai, Xj) : < Xj < Xj-i < ... < Ai < 1} . 



Tr{e~ tHv ~e- tIi2 ) 



pi 2) (o) 



1=1 



J m V(0)d6, ci(V) = J2 F )> C^ = (2n) 



j+n=l 
j>2 



(2 



/ / { Lf (A, 9) \ V{- V^TT ViO^MidXidXj , and 



Lf{x,e) = Y,^- x ^) 



k=l 



i=l 



3-1 fe 

^(A fc -A fc+ i)^i 

fc=l i=l 



In particular, for J = 2, the formula gives 



(1.3) 



Tr(e 



-t-Hv 



p! 2) (o) 



+ t 



y(6»)d6» 



2! 



y 2 (e)de = o(t 3 ), 



as t 4 which is the van den Berg [25] results under our assumption on V. For J = 3, the formula 
gives 



Tr(e" 



-tHy 



-tH 2 



p! 2) (0) 



+ / / v(9)de~ t - I v 2 ((->),w 



(1.4) 



^ / V 3 (9)d9 + ^- I \VV{0)\ 2 d6 = O{t% 



as t I 0. 

For ei = 1, a recurrent formula for the general coefficients in the expansion was obtained in the 
seminal paper by McKean-Moerbeke [19] using KdV methods. Using these techniques, and the 
symmetry of certain integrals, Colin De Verdiere [11] computed the first four coefficients in K 3 . The 
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results in this paper are motivated by [6] where (1.2) is proved by Fourier transform methods for 
all d > 1. Our proof is a combination of probabilistic arguments and Fourier transform techniques 
and unfortunately is much more technical than [6], These results are also motivated by [3] where 
an analogue of van den Berg's results [ 15] (the computation of the first two terms) is proved for the 
fractional Laplacian and other related non-local operators. It is interesting to observe here that 
(integration by parts) 

/ \VV{9)\ 2 d0 = [ -AV(6)V{d)dd = S{V,V), 

JS, d JM d 

which is the Dirichlet form of V with respect to the Laplacian. Based on this, it is natural 
to conjecture that the third term in the expansion for the fractional Laplacian should involve the 
Dirichlet form of V for the operator (— A) Q / 2 . But this is not the case, as we shall see momentarily, 
which is somewhat surprising. 

To state our results for stable processes, we briefly introduce the a/2-subordinators in order 
to more clearly exhibit the similarities and differences from our result to the Banuelos-Sa Barreto 
[6] result. For < a < 2, an a/2-subordinator is an almost surely non-decreasing [0, oo)-valued 
process S = {St} t>0 starting at and uniquely determined by its Laplace transform 

E[e~ xs *] =e~ txa/ \ 

for all t > and A > 0. Throughout this paper we will often write 5i jCt /2 for S± to emphasize 

the a dependence. We write Z = Y for two random variables Z, Y with values in K d to mean 
that they are equal in distribution or have the same law. That is, for any Borel set B C M. d , 
P(Z £ B) = PiY £ B), where Z, Y could be defined on different probability spaces. Our analogue 
result to (1.2) for the fractional Laplacian (—A) ^ is provided by the following theorem. 

Theorem 1.1. Let < a < 2 be given. Suppose V £ S(R d ) and denote the fractional Laplacian 
and its associated fractional Schrodiner operator by H a = (—A) 2 and Hy — H a + V, respectively. 

Denote the heat kernel for (— A)^ byp[ a \x) (see (2.3) below). Assume that M > 1 is an integer 
satisfying M < . 

(a) Given J > 2, there exists a bounded function i?j Q ^ 1 (t), < t < 1, such that 

(1.5) rr(e ~ t t ) ~ e ~ tgn) = -* / v{6)de + y: Y,(-ir +j C%(V)t^ + t»SS.W4(t) 

Pt (0) J* d j=2 n=0 

where 

$^ 1 (M)=rmn|j + l,2 + ^|, 
and the constants C$°j(V) are given by 

<]m = / ( l_ m * [s#> ">} ] n- g « n 

4'>»(A,«) = gs;,_ Al 

fe=l 



Si, 



T d/2 



fe=l 



P [ a \oy 
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where the X' k s are as in (1.1). Moreover, the random variables S Xi _ X2 , S X2 _ X3 ,...,S X ._ 1 _ X . I 
S*_, Xi _ x ,s are independent and satisfy 

j'-l 
fe=l 

and S* = Si for any I e {1 - (A x - Aj), A fc - A fe+ i}^. 
(b) -For any j > 2 and 1 < n < M , 

^ C n}{V) = C^{V). 

We note that when a — 2 the last Theorem remains true and S\ k -\ h+1 = Xk — Xk+i and the 
condition on d and M is not needed. The reason for this is that in the later case, S t = t and then 
Si i = 1. What part (b) in the theorem proves is that our results are robust. 

To see the connection to the Bahuelos and Sa Barreto result more clearly, we state the following 
theorem which is an immediate consequence of Theorem 1.1 (see §8) and which resembles (1.2) 
more closely. 



Theorem 1.2. Under the same conditions of Theorem 1.1, we have 
Tr{er tHv - e- tH ") 



Tr (p- ta v _ P ~iti a \ r . . , , , 

(1.6) lr(e 6 1 = t / V(6)d6 + J2 (-l) n+J C% (V)t^ 

— I /(Hi /rod 

^+j<*S+i(A f ) 

2<3<J, 0<n<M-l 

+ o(i*S+ ) i(^)), 

as 1 1 0. 

Now, to obtain (1.2) from the last theorem we note again that for a = 2 we have no restrictions 

on J and M other than J > 2 and M > 1. Also observe that $ j+ x (M ) = min{J + 1,M+ 2}. 

Then, by taking M = J — 1 wc conclude $yj x ( J- 1) = J+l. As a consequence of (1.6), we arrive 
at 



- = ~t V{6)d6 + V {-l) n+ 3C%]{V)t n+1 + 0(t J+1 ), 

jRd n+ 3 <J+l 

2<3<J, 0<n<J-2 



as t X 0. But, notice that in this case, 



£ (-irM!^)* n+i = E c ^< 

n+j<J+l Z=2 
2<3<J,0<n< J-2 

and (1.2) follows. 

In §11 we provide more specific expansion formulas for a 's of the form 2/k, where k positive 
integer. These examples are the only cases where ^ + j are integers for all n,j, because for the 
particular case n = 1 and j = 2 there exists an integer too > 3 such that — + 2 = too , which implies 
that a = — — o. 

mo — 2 

The assumption M < in our theorem is sufficient to prove two crucial facts needed in our 
expansion. Namely, (1) that the coefficients in Theorem 1.1 are finite and (2) that the remainders 
that appear in the definition of R^^t) (see §4 below) are bounded for t E (0, 1). Since M > 1, 
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the condition M < determines, for a given d, the range of a's for which Theorem 1.1 holds. 
Thus, for example when M = 1 and d = 1, Theorem 1.1 only permits the range 1 < a < 2. In §7 
we will show how a modified version of this condition (namely 4^- — | < can widen the range 
of a's for which Theorem (1.1) remains true when d = 1, 2, 3 and M = 1, 2. 

A particular case of Theorem 1.1 and our results in §7 is the following corollary which extends 
the results in [>] where the second coefficient is computed. 

Corollary 1.1. 

(i) Ford=l , 

Tr ^ tH l~f tHa) + t f V{6)d6 -trJ v 2 (9)d9 + £ [ V*(9)d9 
Pt (0) J R 



0(t 2 +^), a/a 6 (1,2), 
0(t 4 ), if a G (1/2,1], 



as 1 1 0. 

(ii) For d — 1 and | < a < 2, we have 
Tr(e- tHv - e- tH <*) 



as 1 1 . 



Pi '(0) 



+| y F 3 (#)d# + C ha t 2+i J \W(9)\ 2 d9 = 0(t 4 ), 



(iii) For d > 2, 



] +t V(9)d6-- V 2 (9)d9 + - I V 6 {9)d9 



0(t 2 +i), if a 6(1,2), 
0(i 4 ), a 6(0,1], 



as t J, 0. 

(iv) For d>2 and 1 < a < 2, 



.Ld w 2! 



+ S / 7 3 (0)d0 + £ d , Q f 2+ ^ / |VV^)| 2 ^ = 0(t 4 



as t J, 0. 

(v) For d > 3, | < a < 1, 



^ a) (o) 



+ tf V(e)d9- t - f V 2 (9)d9+ t - f V 3 (9)d9 

JR d z - JR d <J ; JR d 

^»d0 + £ d , Q i 2+ * / \VV{9)\ 2 d9 = 0{t b ), 



4! 



as £ \. 0. 
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Also for d > 3 and \ < a < 



Tr(e 



-tHi 



+ t [ V(9)d8- t - [ V 2 (9)d9+ t - I V 3 (9)d9 



as t 1 . 

(vi) For d > 4 and < a < \, , 



2! . 

- I V i (9)d9 = (D(t 5 ), 



3! 



Tr(e 



-tH v 



-,-tH a 



P { t a \o) 



+ t 



V(9)d9^ t - f V\9)d9+ t - 
^£v*(9)d9 = (D(t 5 ), 



V 3 {9)d9 



as t I . 



The constants £d, a are defined as follows: 



_ C d ,gKi(d : a) 



T d/2 



(27T)° 



p[ a \oy 



with 



KM a) 



l r \ 



"'O 



E 



^1— w^w 



dwd\\. 



The question of whether our result holds regardless of the choice of d and M as in (1.2) remains 
an interesting open problem which reduces to verifying that the expectations in the formula for 
C^*j (V) are finite for all n and d. 

To gain a better understanding of the applications of the robustness result (part (6)) in Theorem 
1.1, which is proved by means of weak convergence, consider the following special case of Corollary 
1.1. For all d > 1 and | < a < 2, we have 



Tr(e 



-tH v 



-tH a 



-+t 



V(9)d9 - - 



V 2 (9)d9 



4 / V 3 (9)d9 + £ d ^+i [ |W(0)| 5 



d9 = 0(t 4 



as 1 1 0. Interestingly, due to part (b) we see that £d,a -> ^ as a | 2, despite of the fact that thus 
far we are only able to provide a representation which enables us to conclude that £d. a arc hnitc 
and strictly positive with no other explicit knowledge for this quantity. 

The rest of this paper is organized as follows. In §2, we introduce notation, definitions and 
probabilistic facts about stable processes and subordinators. In §3, we find formulas (and bounds) 

for the difference of the Fourier transform of the heat kernels p^ v (£, rj) — (£, rj) that will allow 
us to express the trace in terms of Fourier transforms. In §4, we prove the boundedness of the 
(J+l)-th term in the trace formula found in §3 for t £ (0, 1). In §5, we simplify the trace formula 
given in §3 by finding an explicit value of certain integrals in terms of the subordinator process. 
In §6, we use an elementary Taylor expansion formula for the exponential function to define the 
coefficients and remainders involved in Theorem 1.1. In §7, we give an improvement of Theorem 1.1 
when d = 1,2, 3. The proof of Theorem 1.1 is given in §8. In §9, we compute some coefficients and 
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provide a representation for the constants Ld, a - In §10, we prove Corollary 1.1. Explicit expansion 
and examples are provided for some particular a's in §11. Finally, in §12 and §13, we extend 
the result i) and hi) in Corollary 1.1 to a-relativistic processes and to mixed-stable processes, 
respectively. 

2. Stable Processes and subordinator. 

Let X = {X t } t>0 be the d-dimensional symmetric a-stable process of order a € (0,2]. The 
process X has stationary independent increments, which means that if s, t > 0, then the increment 
X t + S — X t is independent of the process (X u , < v < t) and has the same law as X s . Moreover, its 



transition density Pj {x, y) = pf^ (x - 
(characteristic function) as 



y), t > 0, x, y € K d , is determined by its Fourier transform 



-t|?l c 



E°[e- 



i<£,X t >] 



-i<jl,f > ( Q ) 



Pi \y)dy, 



for all f > 0, £ G M. d , where P x and E x denote the probability and expectation, respectively, of the 
process starting at x. Then, we have for any Borel set A C K d , 



P x (X t eA)= f p i f\x~y)dy. 

J A 



Henceforth, E will denote the expectation of both an arbitrary random variable or processes started 
at 0, whereas a.s will mean almost surely. 

We recall again that an a/2-subordinator is a.s non-decreasing [0, oo)-valued process S = 
{5f} (>0 which also has stationary, independent increments, starting at and uniquely determined 
by its Laplace Transform 



E [e 



-IX 



,/2 



for all i > and A > 0. 

Notice that the last equality implies that for all £ € M. d , 



(2.1) 



E 



-l?l 2 s t 



-*l?l c 



We point out that this equality is the link between the results in this paper and the results in [f 
and it will be used several times throughout the paper. 

It is a standard fact (see [8, p. 22]) that the a-stable process X can be obtained as a random 
time change of Brownian motion where this random time is an a/2-subordinator. In other words, 
we can write Xt = i?2S t where B is a d-dimensional Brownian motion and S is a a/2-subordinator 
and these are independent processes. Using this last fact, the Fourier inversion formula (see §3), 
we see that 



(2.2) 



p[ a) (x) = (2ny 



e i<£,*> e 



= E 



-d/2 f 



where p[ 2 \x) = (4irt) d / 2 e « and ^"'^'(s) correspond to the transition densities of the d- 
dimensional Brownian Motion process {i?2t}t>o and the a/2-subordinator, respectively. It follows 
from (2.2) that p^ (x) is radial, symmetric and decreasing in x. Moreover, these functions satisfy 
the following scaling property and inequality. 
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(2.3) 
where 



(a), n , w d T{d/a) 
Pi (°) _ 



(2ir) d a 

Here and for the rest of the paper, Wd denotes the surface area of the unit sphere in R d . From 
equation (2.2), we also conclude that 



(2.4) (4ir) d / 2 p[ a) (0) = E\s-y 2 

In fact, we claim that for all — oo < r\ < j, 

r(i - 22) 



(2.5) 



E 

a/2 



r(i - rj) 



To see this, we observe that ( ZS 1 a) = Z, where Z — exp(l) is an exponential random variable 
with parameter 1 independent of S^s, . By independence we have 



P([ZS-i) <A) = 



e- u du } ri[ a/2) (s)ds 



1 - E 



1 - e A = P(Z < A). 



Hence, it also follows by independence that 



E [Z-v] E 



= E 



provided E {Z~^\ and E 



are both finite. But, this only holds when — oo < i] < 2., since 

/•oo 

< E [Z' 1 ] = / *(T +1 )- 1 e-'da = T( 7 + 1) < oo, 
Jo 

when 7 + 1 > 0. 

The equalities (2.4) and (2.5) will be useful in proving the finiteness for the coefficients and the 
boundedness of the function i?j Q ) : 1 (t) in Theorem 1.1. We also mention that the condition given 
in Theorem 1.1 is derived from (2.5). 

3. Heat trace in terms of Fourier transform. 
Let V denote the Fourier transform of V £ S(R d ) with the normalization 



v(0 



-i<x,{> 



V{x)dx. 



We note that because of our definition of V, we have 
(i) (Inversion formula) 



V(x) 



1 



(27r) d 
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(ii) For f,geS(R d ), 
(3.1) [ e -i<*,e>f( x )g( x ) dx= l I jWWi.-<ruW. 



(27r) d 

We recall that the linear operator H a ~ (—A)"/ 2 with domain 

{feL 2 (R d ) : \£\ a f(OeL 2 (R d )}, 

a e (0,2), is defined via Fourier transform by H a f(£) = |£| Q /(<!;)• This operator, often referred 
to as the fractional Laplacian, is essentially a self-adjoint operator on Cfi°(R d ) with spectrum 
Spec(H a ) = [0,+co). For V € S(R d ), we call H v = (-A)f + V the fractional Schrodinger 
operator with potential V, where V acts as a multiplication operator. Hv is self-adjoint in the 
domain of H a with Spec(Hy) = [0, +oo). We write e~ tHa and e~ tHv for the heat semigroups 
acting on L 2 (R d ) with heat kernels given by p\ a \x, y) and p? v (x,y), respectively. We refer the 
reader to [24, 13] for general definition and spectral properties of the semigroup e~ tA of the self- 
adjoint operator A on Hilbert-space. For our purpose in this paper we recall the Feynman-Kac 
formula which gives the heat kernel for Hy ■ That is, 

(3.2) p^{x,y)=p { t ) {x,y)Ei. 



e -f*V(X s )d s 



where E\. is the expectation with respect to the stable process (bridge) starting at x and condi- 
tioned to be at y at time t. For more detail about formula (3.2), see [3, 24]. 

Our goal now is to derive a formula for Tr(e~ v — e~ tHa ) for the fractional Laplacian similar 
to the one in [6] for the Laplacian. 

Proposition 3.1. Let V € S(R d ), then 

T ^-' Hv --*> - (41 (^<?.-«>-^(?.-«>) * 

Proof. For alH > and x, y £ R d , we have 

(3.3) d tP { t a \x,y)=-(-A)!p{ a \x,y) 
and 

(3.4) d t p? v (x, y) = - [(- A)| + V{x)]pf v (x, y) 

Po v (x,y) = S{x - y). 
By taking Fourier transform on R 2d , we deduce that 

(3.5) pl? ) (Z,v) = W 1 S(Z + ri) 
and that 



(3.6) P£ V (Z,V) = (27r) d 6( V + 0. 

Now, by directly solving (3.5) and (3.6) we find that 

^ ) (^r 1 )^(27T) d S(f 1 + Oe- mC 
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and that 

(3.7) pP«,»?)-pP&»7) = - i /* f e-^-^ a v(e)^^-e, v )ded S . 

(27r) d J J Rd 

On the other hand, from (3.3), (3.4) and Duhamel's Principle we see that 

p? v {x,y)-pf\x,y) = - [ f p[ a J s {x,z)p^{z,y)V{z)dzds. 

J0 JR d 

Hence 

Tr(e-* H v - e- tH ") = / (p? v (x,x) -p[ a \x,x)) dx 

JM d V 7 

(3.8) = - I I Pt-s(x, z)Ps V {z,x)V(z)dzdxds. 

J0 JR 2d 

Expressing the right hand side of (3.8) in terms of Fourier transform we arrive at 

Tr(e- tHv - e- tH ") = -—^ [ f V{9)plF 'faMl^-p - 9, -r)d^drd9ds. 

Since 

P&(-H-6,-t) = {^) d 8{r + 9 + fi)e-^-^ e+ ^ a , 

we see that 

Tr(e~ tHv -e- m ") = --^ I I e^-Wrf* (-r - 0,T)V(O)dTdOds. 

( 27r ) Jo Ju 2d 

The conclusion of the proposition follows by setting £ = — r in last equation, r/ = £ in (3.7) and 
integrating with respect to £. □ 

If we now iterate the equation (3.7) J-times, we obtain 
Corollary 3.1. Let V G S{R d ), < a < 2 and set 

where s — (s±, Sj), Sj, < Sfc+i and = (#1, 0j-i). Then for J > 2, 

(27r) d y 7 Rd 



T /"••• f ' / F^\s,^ P if^-j20i,v)f[V(9 i )d9 i ds i + 

3=2 { ' J ° Ja Ja jRjd i=l i=l 

/•* /-si r SJ r 1 \ J+1 ,7+1 ^ 

II- / ^\( s ,e,%f/ +1 (C-E^^)Il^) d ^ 

Jo J'R( J + 1 ) d ^ f = l 



( 27r )(J+i)«« 7 7 
Furthermore, we conclude 
(3.9) Tr(e- tHv - e~ m ") = -^ a) (0)F(0) 



j=2 V ' JO JO JR3 d i=1 i=1 

(-t) /+1 z" 1 /" Al /" Aj /" -= — - — 1 — 1 



(27r) (.7+2)d ^ ^ 



Jo Jr(-'+2)« j+1 f-f M 
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4. BOUNDEDNESS OF THE (J + 1)-TH TERM. 

Our goal in this section is to provide an upper bound for the absolute value of last expression 
in (3.9) in terms of p[ a \o). The following Lemma is a consequence of (3.1), the inversion formula 
and induction. Therefore its proof is omitted. 

j J 
Lemma 4.1. Let J > 1 and {pi} i=0 C <S(R d ) radial functions. Set = Po{0 II Pjilj — 0> 

3=1 

where jj € M d are constant vectors. Then, 

^ p -i £ (7jiij) 3 3 

*j(7) = (27r) d / e '= l " ' paC^Xj - 7 ) \\ P] (x 3 )dx 3 . 

Remark 4.1. // we set po — p[?i_ x \ and pj — p|(a--a- ) * n ^ e ^ as ^ Lemma, with jj = ®k> 

J fe=i 

we obtain 



(4.1) 



(2tt)° 



i=i j'=i 

Next, it is known that the transition density p^ (x, y) satisfies the Chapman-Kolmogorov equa- 
tion, namely, 

(4.2) / p(f\a-z)pi a \z)dz=p[t(a), 

for all a £ M. d and t,s > 0. With this equality at hand, it easily follows that 
r J J 

( 4 - 3 ) / piu-AoE^-^-^nptri-Aj+o^j)^ =Pt(Lx J+1 )( x -y)- 

jRjd 3 = 1 3=1 

It can also be proved by means of the inversion formula that 

r — 1\ («,7J+i> + E (1j> x j) \ J+1 3 3 

(4.4) / el ^ ' J ^y(w = (2^) ( ' /+1) M^)^^(-I>^ + 2; )• 
• / K <J+l,d f=i r=i ^ 

Proposition 4.1. Assume < i < 1 and define 

rj+i(t)= f t 1 ■■ r [ F^(t\^,e) - 'j2o i ,-OY[v(e l )dmd\ i . 

Jo Jo Jo Jm<J+2)d + ~[ fJl 

There exists a positive constant C = Cj+i d a (V) such that 

\r J+1 (i)\<Cp ( t a \0). 

3 

Proof. Set 7j = #i and 

J J 

P (x - v,{*3}i=i)=pi(i-\ 1 )(E x 3 -( x ~ y)) n^ri-A J+1 )(^)- 
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Then, by definition of Fourier transform and (4.1), we have that 
I J+1 (tX, 6) = [ F(% (tX, 8)p^ i (£ - 7./+1 , ~m 



F^%{tX,^6) / e -m,*-y)+(*,-tJ + i)} p fv f Xty)dxdyd £ 
p^ j+i (x,y)e-^-^^F { J %(tX,x- y,9)dxdy 



=(2?r) e 1 ) P t \ J+ i( x ^v)p{x-y,{x j }- 1 )dxdydx ] . 



Jd J R 2d 



Now, because of (4.4), (4.3) and the fact that pf£ i (x, y) < e t||y|l ^°°(K d )p[^ i (x, y) (which follows 
from (3.2)), we obtain from the last equality that 



J+i 



I J+1 {tX,9)T[V(6 i )d6 i 

J 

(2^+2)^ / JJ V( _ x, - x)p$ (x, y)p(x - y, {x } J )dxdydx 3 

(2 7 r)^ +2 ) d ||l/||^ (R£i)e t||y|1 -^) / \V(-x)\ [ p { «l_ x Jx-y)p^(x-y)dyd, 



< 



It follows from (4.2) that 

|r w ml ^ ) ^f ( ° ) IWIi-^ l ""'-*-llMI W 

As a consequence, we have that 

rj+i(*) 



□ 



i2j+i(t) = 



(27r)( J + 2 )^ a) (0) 
is bounded for < t < 1. 

5. Heat trace computation by means of subordinators. 

In this section we further investigate formula (3.9) involving TV(e~ v — e~ tHa ). We start by 
integrating the function 

(5.1) Fj Q, (tA,^0)e -u ^ fl =e 

with respect to £ over R d . The integral could be easily computed when a — 2 using two elementary 
facts. Namely, for any 7 G E d , 

(5.2) |e + 7| 2 H£| 2 +2<£,7>+l7| 2 
and 

(5.3) / e-^-^dC = ^'H- d ' 2 . 
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Unfortunately, we cannot calculate (5.1) in the same way because there is not a close form for 
|£ + 7| Q , when < a < 2. Instead, we will follow a probabilistic approach by means of a/2— 
subordinators and their Laplace transform given in §2 that relates | • \ a to | • | 2 to find the value of 
the integral involving the quantity in (5.1). We begin by observing that (2.1) implies that for all 
c > and t > 0, 



(5.4) 



In addition, for any sequence of numbers {Afe}^. =1 , j > 2, satisfying 
(5.5) < Xj < Xj-i < ... < A 2 < Ai < 1, 

we have 

3-1 

&i = i ~'(i—{x 1 —Xj})+\ 1 —Xj ~ Sx^-xj + ^2 (S Xk -x k+1 +(x k+1 -x 1 ) - Sx k+1 -X;j) 



fc=i 



For 1 < k < j ' — 1 consider the random variables 
and 



S*x k -x k+1 — Sxk-Vn+CVu-A^) Sx k+1 -x 3 



■ST-(Ai-Aj) - "^l-fAi-A^+Ai-Aj _ S'Ai- 



A,- 



Since the process S has independent and stationary increments, we see that the random variables 
{^Afc-Afe+i) "^*-(Ai-A )} are independent and furthermore, 



S 



A fc -A fc 



= <7 



(5.6) 

We also have, of course, that 
(5.7) 



1— (Ai— Aj) - ^l-(Ai-Aj)- 



3-1 



S, *-(A 1 -A J ) + ^2 S **k-X k+1 - Si- 



fc=l 



As before let us denote, for simplicity, jk — ^ follows from (5.4), (5.6), (5.7) and the 

i=l 

I 3-1 



independence of l^k-Ajs+i' "^*-(Ai-A )} that 

3-1 

e _ t (i_[ Al -A J ])icr Y[ e -*( A *-^+i)ie-7*i 
fe=i 



(5.8) 



E 



E 



E 



exp -* 2/a E^-A fe+ je 



Tfcl 



fe=i 



(^ 2/Q j^-CA.-A,)^! 2 +E 5 A fc -A fc+1 l^7 fc | 2 } 



Next, consider the random variable 
(5.9) 



Lf(X,0)=Y,Sl 



Afc — A fc + 1 l7fc| 



fc=l 
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where A = (Ai, Xj) satisfies (5.5). By (5.2), (5.7) and completing squares we easily get that 



5r_ (Al _ A3 .)iei 2 +E^-A fc+ J^-^i 2 = 5 a 



fc=l 



1 ^ 

fc=l 



+ 4 Q) (A,0). 



Also, observe that by (5.2) and the scaling property (2.3) we have 
exp | -t 2 /"5i 



1 J " 1 



fc=l 



(5.10) 

over the the set where < Si < oo. Here, 
(5.11) C djl 



= c d ,^ tt) (o)5r d/2 , 



r d/2 



We now combine these calculations to find the value of the desired integral. More precisely, we 
have 

k 

Lemma 5.1. Let A = (Ax, Xj) satisfy (5.5), — X) @i an d ® = (^i> •••) @j-i)- Then, 

-(a) 



L^(A,fl)>0, a.s. 



(5.12) 



-d/2 - t 2 ' a L ( f ) (\,e') 



Proof. Assume jtc = (bi,k> ••-> &<2,fc)- By Cauchy-Schwarz inequality and (5.7) 



3-1 

fc=i 



= \j2 S k-X k + 1 b mM> = \j2{ S *X k -X k + 1 } W{ 5 A fc -A fc + 1 } > 

rn=l {k=l ) m=l \k=l ) 

d fj-1 } 3-1 3-1 

m=l Lfc=l J fe=l fc=l 



In fact, under the convention that ^ = 0, we have for j > 2, 



(5.13) L { f } (X,9) = S^ 1 



r=l 
3-1 



3-2 3-1 



^-(Ai-A^E^A.-A.+ijTfcl 2 +E H S, A r -X r+I ^A.-A <>+ xl7r 



fe=l 



r—1 s—r^l 



On the other hand, (5.12) follows by integrating (5.8) with respect to £, applying Fubini's Theorem 
to (5.8) and using (5.10). □ 

Remark 5.1. We note that from (2.2) and the fact that < S\ < oo, a.s, 



(5.14) 



< E 



q -d/2 .-t'/UflV) 



< £ 



5' 



-d/2 
l,a/2 



< OO. 



In fact, all the above results are true for a = 2 in which case Si t i = 1 and all our calculation 
considerably simplify. 
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6. Bounds for remainders and coefficients. 

We observe that the exponential function is involved in (5.12) and that this term is part of 
the expression for Tr(e~ v — e~ tHc ") in (3.9). Our next step is to use a Taylor expansion of the 
exponential function with a particular remainder to obtain a finer estimate for the trace. This 
implies, as the reader may note, that in (5.12) we will have to deal with expectations. Hence our 
goal in this section is to give conditions to guarantee the finiteness of these expectations. Once 
this is done, it will follow easily that the coefficients and remainders to appear in (3.9) are also 
finite and bounded, respectively. 

We recall the well known expansion for the exponential function 



(6.1) 



y(-i)y , (-1)' 

n=0 



valid for every x > and m > 1, where we call (3 m (x) € (0, 1) the remainder of order m. With this 
expansion at hand, we now introduce two functions from which we will obtain the desired finer 
estimates in the trace formula (3.9). For j > 2, 



(6.2) 



Jo Jo Jo Jmi d i=1 i=1 



c d , aP { t a) (o) 

lo Jo Jo 
The remainder function is 



1 /-Ai 



E 



a-d/2 -taiW(A,0) 
,5 l,a/2 e 3 



J-l J-l 

i=l i=l 



Cd,, 



ml(2ny d J Jo 



1 



Xj-i 



(6.3) 



:V(-J2^)Y[V(^)d0 t d\ t d\ J , 

i=l i=l 



where the random functions f3* mJ (t) = t 2 / a L ( *\\9)[3 m {-t 2 / a L ( f' \\,6)) are nonnegative. 
Remark 6.1. We note by (5.14) an d (6-2) that 

-d/2~ 



T d {j,t) 



Cd, a E 



< 



S 



l,a/2 



J-l J-l 

n-E^)IP(^) 

i=l i=l 



j!(2vr)^ J R(j - 1)d 

for all j > 2, d > 1. Observe that the left hand side is finite since V G 5(R d ), proving at the same 
time the finiteness ofTd(j,t), for all t > 0. 



We now proceed to prove the finiteness of the remainder-functions in (6.3) and define the 
coefficients given in Theorem 1.1. 

Lemma 6.1. Assume M > 1 is an integer satisfying M < . Then, for all t > and j > 2, 



R„d (*) ^ C 3,d,M 



El 



M 



7k I 



j'-i j-l 

n-E^)IP(^ 

i=l i=l 



Mi, 
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whe 



C 



Cd.aE 



j,d,M 



s 



M-d/2 
l.a/2 



j!M!(27r)J' d 



i=l 



Furthermore, 

(a) for all integers n < M, 



(b) and for all j > 2, 
T d (j,t) 



0<E 



M-l 



?-f {Lf\\,e)}' 



< oo, 



j2(-irci:](v)t-+(-irt^R%\t), 



where 



c ( n a kv) = 



1 



Xj—i 



E 



5 W2 {4 Q) ( A ^)}"] II V(^)dO l d\d\ 



i=l i=l 



< OO 



Proof. We start by observing that the condition M < guarantees that < E Jsj 1 /' 2 

for all integers n < M, according to (2.5). From this we proceed to prove (a) as follows. Recall 
that 

i-i 
fe=i 

It follows from (5.9) and last equality that 

< Lf{\,6) < XX_ Afc+ >| 2 < S ha/2 J2 \j k \ 2 

k=l fc=l 

Then, from the last inequality we conclude that for all integer n < M, 

'i-i \ 71 



E 



S 



< E 



S 



t-d/2 
l.a/2 



\k=l 



Thus (a) now follows easily from last inequality. On the other hand, (b) follows from the Taylor 



expansion (6.1) applied to (6.2). We remark that C ( n a Uv) = i?|"J(0). Therefore the last expression 



is finite, according to (6.1) 



□ 



Remark 6.2. Lemma 6.1 allows us to bound the remainders by a constant for all t > and shows 
the finiteness of both the remainders and the coefficients C 1 ^"- (V) by proving the finiteness of the 
expectations under the condition n < M < ^rp- . Indeed, this condition is introduced to make sense 
of the Taylor expansion of order M when it is applied to the function (5.14). As a consequence 
our results are dimensional dependent. The reason why this does not happen when a = 2 is that 
in this case the time change is trivial, St t \ = t, and L^p is nonrandom function. These two facts 
considerably reduce all above computations, thereby the dimension only appearing in the integrals 
involving V, which are finite since V £ S(M. d ). 
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7. AN IMPROVEMENT FOR DIMENSION d — 1,2,3. 

We recall the basic inequality 



O'-l 



M 



3-1 



EH <o--i) M "E«^ 



fe=i 



,-Af-4 



7* 



1 2 M 



valid for all j > 2 and positive numbers {cifc}{ =1 ■ 

From last inequality and (5.13), it follows for all integer M > 1 that there exists a constant 
Cj,M > such that 

L - 1 fe=l L 

(7-i) E E ^[(5t-A,. +1 ^ s -A s+1 ) M ^r 

r—l s— r+1 

whenever the expectations involved in the last expression arc hnitc. The purpose of this section is 
to provide conditions under which these last expectations are finite for dimension d = 1, 2 and 3. 
We proved in §5 that 

j-i 

•ST-CAi-Aj) + E ^Afc-Afc+i = 



l7r - 7s I 



2jU 



fc=l 



for {Afc} fe=1 satisfying (5.5) and where the random variables on the right hand side of the last 
equality are independent. In particular, it follows that Si > S£ + Sf for any distinct lo,h G 
{1 — (Ai — Aj),Afe — Afc+i}^!- Now, observe that each expectation in (7.1) can be written as 



E 



(sis; o ) M s; 



M- 



and these expectations satisfy 



E 



{s* h si 



* \M 



Si 



M+l 



< E 



* Q* \M 



Lemma 7.1. Let j > 2 and {Afc}^, =1 satisfying (5.5). Let Iq,Ii be two distinct numbers 



{1 - (Ai - Xj), A/s - Afe + i}£ =1 . T/ien, 



< E 



< 00, 



provided that M/2 - d/i < a/2. 
In particular, when 

i) M=l, if^ <a. 

ii) M=2, if±4<a. 

Proof. Because of the inequality 2{ab) 1 / 2 < a + b, for any a, b > 0, we have that 



M 



{Sto+St JM +d/ 2 



< 2 - M - d ' 2 E 



(st si) 



M 



(Sf Sf) M / 2+d ^ 



Now, recall that S* and 5 ; * are independent and S* = l 2 / a Si a / 2 ■ Therefore 



(7.2) 



E 



(5; o s , *) M /2+<i/4 



. 2 fM dl / 



M d 

S*<* 4 
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The result follows from the inequality M/2 — d/A < a/2 which guarantees the finiteness of the last 
expectation. □ 

As an application of Lemma 7.1, we have 

Corollary 7.1. Assume j > 2. 

(i) For i < a < 2 and d — 1, we have 



E 



S-i /2 £f(A, 



< oo. 



(ii) For d — 1 and | < a < 2, d — 2 and 1 < a < 2, d — 3 and j < a < 2, we have 



E 



S 



< oo. 



The following is a version of Lemma 6.1 for dimension 1,2, and 3 where the condition M < ^±2- 
is replaced by -1 < M/2 - d/A < a/2. 

Lemma 7.2. 

(i) For d—1, and M=l, we have {or all | < a < 2 and j > 2 that 

Ti(j,t) 



(2^p| a) (0) 



(ii) For M=2 and j > 2, we have 
T d (j,t) 



(2tt)^ q) (0) n=0 
when d—1 and | < a < 2, d — 2 and 1 < a < 2, or d — 3 and ^ < a < 2. 



where the remainders 



are bounded by a constant for all t > 0, according to Corollary 7.1 
and the fact that M/2 — d/A > — 1 for all M, d as stated above. 

Proof. We start by recalling that for any j > 2 

Ij = {A = (Ai, Xj) : < A j < < ... < Ai < 1} . 



Next, under the notation given in Lemma 6.1, we have that 
are bounded by 



C[ a hv) and R\ M d \t) 



M = 1,2 



E 



s 



- d/2 {Ll a \x,e)} 



M 



i=l i=l 



dXjdXidOi. 



This last expression is also bounded, based on the facts given at the beginning of this section, up 
to some positive constant by terms of the form 



I {l h)^-^dXjdXi (e 



Si\ 4 

1 1 o 



3-1 

Ei 

k=l 



,2M 



J-l 3-1 

II ^ 

i=l i=l 



Mi 



where lo — Iq(X) and h = h{X) are two distinct numbers in {1 — (Ai — Xj), Xk — Xk+i} J k _ v 
Now the term 



(loh^-iUXjdXi = 



1 ,-A, ,.A,_ 



JO JO 



(l l 1 )i(%-i)dX j ...dX 1 
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is clearly finite when M/2 — d/4 > 0, which is the case for M = 1, 2, d = 1, 2 and M = 2, d = 3. 
But, when M = 1 and d — 3, we obtain —1 < M/2 — d/4 = 1/2 — 3/4 = —1/4 and this case deserve 
special attention. 

We observe that for all 1 > \ > Xk+i we have 

Afc+i(A 4 - A fe+ i)~~<£\ fe+ i < / (Aj - X k +iy^dX k+1 < / w~^dw = - - 

7n ,/n 2a -1 



, . , i . , i _ i 2a 



and 

,' , " I (l-A 1 +A J )-^(A J -_ 1 -A J )-*dA i <(l-A 1 )-^ 
Jo 2a - 1 

provided that a> \. Then, it is not difficult to see that 

2a 

(loli)~^d\jd\. t < 

h 



2a- 1 



□ 



8. Proof of Theorem 1.1. 
Proof of part (a): Recall that, for J > 2, we have defined 

R J+ x(t) = rj+l(t) f , 

and also showed, according to Proposition (4.1), that this remainder is bounded by a constant for 
< t < 1. Also M < implies that ii^eP (*)>•••? ^(m+1) anc ^ R-M+2{t) are, according to 
Lemma 6.1, bounded by a constant for < t. 

Next, (a) follows by substituting the terms found in Lemma 6.1 into (3.9). More precisely, 



Tr t-tH v _ e -tH al \ r J^ A S^ 



^ a) (o) 



j=2 n=0 

(8.1) (-Q'+^+xCt) + J2(-l) j+M « + ™R%\t). 

Therefore, it suffices to take 

flW (t) =r *!S»w |(- i )^ii? J+lW + ^(-i)^M^ +2 M i? (M) w | | 

and this proves Theorem 1.1. 

Moreover, Theorem 1.2 now follows by noticing that 

J M-l 

£ £ (-ir+^g(F)t^ = £ (-i) n+j 'cS(v)t^ 

i=2 n=0 M+j^^W^M) 

2<3<J, 0<n<M-l 

+ (-i) n+J 'cg(v)t^+J 

2<j<J, 0<n<M-l 

and 

{T+i = 0(t $ ?iW), 
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as 1 1 0, provided that ^ + j > $^(M). 

Proof of part (b): We begin by recalling several basic facts about weak convergence. A 
convenient reference for this material is [7]. Let X r , X be /c-dimensional random vectors, possibly 
defined on different probability spaces. We recall that X r converges weakly to X, denoted by 
X r =*> X, if 

lim F Xr {x) = F x (x), 

r— >oo 

for every continuity point x of Fx, where Fx r and Fx are the distribution functions of X r and X , 
respectively. The followings statements are consequences of weak convergence and we state them 
here as a facts, referring the reader again to [7]. 

Fl Suppose that h: R k —> W is measurable and that the set Dh of discontinuities of h is 
measurable. If X r X and P(X e D h ) = 0, then h(X r ) =>• h{X). 

F2 Let X r ,X be A:-dimensional random-vectors. Then, X r =£> X if and only if for every 
bounded, continuous function /, we have 

lim E[f(X r )]=E[f(X)}. 

r— >oo 

F3 Let Y r ,Y be real valued random variables. If Y r =>■ Y and the {Y r } r&N ^ s uniformly 
integrable, then Y is integrable and 

lim E [Y r ] = E [Y] . 

r— ¥ oo 

It is a standard fact that if 

su V E[\Y r \ p ] < oo, 

r 

for some p > 1, then {l^} rgN is uniformly integrable. 
F4 Let X r , X be fc-dimensional random vectors. Then, X r =^> X if and only if for all t> € 
< u,X r > < i>,X > . 

Since we are only interested in a's close to 2, it suffices to prove that if n, d > 1 are positive 
integers satisfying n < ii= , then for all j > 2 we have 

lim C<%\v)=C™(V), 

r— too u u 

for any sequence {a r } rgN satisfying | < a r < 2 and a r t 2. 

To prove last statement, we need to introduce some notation. We recall that Ij C MP has been 
defined as 

Ij = {A = (Ai, Xj) :0<Xj< Xj-i < ... < Ai < 1} . 

We also set 

X r (X) = (^_(x,-^),fii i ^j.j-A,,^ i ^Ai-Aj a), 
X(X) = (1 - (Ai - Xj), Xj-i - Xj,..., Ai - A 2 ), 



3-1 3-2 3-1 

^0 J2 Xk\~fk\ + Y2 J2 X m X s \j m -J 3 \ 
k = l m = l s = m+l 

n 



h n ,d(xQ,xi, ...,Xj) = < 



3-1 

12 x k 

k = 



i+4 



-, for Xq > 0,..., Xj > 0, 



otherwise. 
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With this notation, 



K. d {Xr{X)) =^i 2 {4 ar) (A,0)}' 



We now divide our proof into 5 steps. 
Step 1. X r (X) X(X). 

To see this, we recall that for t e {1 — (Ai — Xj), Xk — Afc+i}^* and A > 0, 

tAQr /2 



E 



This expectation corresponds to the Laplace transform of S£ Q y 2 and uniquely determines its 
distribution. We conclude that 



lim E 

r— f +oo 



e - A5, t*Q r /2 



-/A 



Thus, Sj a j2 =>■ i- On the other hand, due to the independence of S'*_^ 1 _ A .) a / 2 > _i— A q /2'- ■ • 
^Ai-A 2 " /2> ^ e ^ ac * = e " * s bounded in R and F2, we also obtain for every v £ W +1 

that 



-i<t>,X r (A)> 



-i<v,X(A)> 



as r — > oo. The result follows from F4. 



Step 2. V*(-MA)) => K A (X(X)). 

We note that each component of the vector X(X) is positive. Thus, by our definition of h n ^ it 
is clear that X(X) belongs to the set of continuity points of h Hl d- Then, P(X(X) E Dh n d ) = and 
the result follows from Fl. 

Step 3. {h n: d(X r )} reN is uniformly integrable. 

We shall show that there exist ap> 1 and a function C(n, d, p, 9) > such that 



.2) supE[{h n4 (X r (X))} p ] < C(n,d,p, 

r 

To do this, we consider two cases to determine a proper p. 

Case 1. Suppose n — | < 0. In §6, we proved that 

d /i-i \ n 



To prove (8.2), it suffices to show that sup E 
32 that 



is bounded for some p > 1 . Recall from 



(8.3) 

provided that 
(8.4) 



< E 



S 



,p{n-i) 



y(i- ^ 



I)) 



T(l-p(n-D) 



< oo, 



p\n~-\<-<l. 
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If n 



0, we take C(n,d,p,t 



EM 



pi; 



and any p > 1, since clearly in this case 



5 



p{n-i) 



If n — 5 < 0, it is clear that (8.4) is satisfied for any p > 1. But, we wish to pick p so that (8.3) 
is uniformly bounded in r. To do this in a suitable manner, we require the following well-known 
property of the the gamma function(see [2.3]). There exists S (1, 2) such that T(z) is decreasing 
on (0,/Lto] and increasing over (/xo,+oo). Now, from the last property and the fact that each a r 

satisfies 1 < < 2, it follows that for any p > max jl, 



(8.5) 

Therefore, 



2p ( d 



r(2) <r i + -ii --n <r i 



2p ( g - « 



sH 1 ,E[{/,„., / (A',(A))f] < Eg 2p } n (]T>| 



pn 



r(i _„(„_!)) 



C(n,d,p,6). 



\k=l 



Case 2. Suppose n — | > 0. Because of the inequality < n — | < |, which is equivalent 
to < 2n — d < 1, we conclude that d = 2n — 1, since (i and n are positive integers. Therefore 
n — | = i. In this case, from the tools we developed in §7 we obtain that < h n .d{X r (X)) is 
bounded, up to some positive constant, by a finite sum containing terms of the form 



(8.6) 



3-1 

fe=i 



r /2^l 1 ,a r /2 



(f-#) 



for any two distinct numbers Zq, ii in {1 — (Ai — Xj), Xk — Xk+iYf,^- 

■2 

Next, due to the fact that S* «,., i = 0, 1 are independent and have law l" r S\ j(Xt /2, < k < 1 
and n — | = |, we conclude that E [/i„ j[ j(X t .(A)) 2 ] is bounded up to some positive constant by 



E 




We also know that 



E 



? i/2 

l,O r /2 



r(i 



On the other hand, the function T{z) is decreasing over (0,1). Next, we observe that each a r 
satisfies 

1 1 1 
- < 1 < -, 

3 ~ a r ~ 2' 

which yields 



For this, 



C(n,d,2,9) = C 




for some C > 0. 
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Step 4. lim E [h n4 {X r {\))} = E [h n , d (X)}. 
This is a consequence of Steps 2, 3 and F3. 

Step 5. Notice that by Holder's inequality and Step 3, we have proved that for some p > 1, 
(8.7) aupE[h n , d (X r (A))} < sup(E[h n , d (X r (\)) p })» < {C(n,d,p,9)}$ , 

r r 
X 

where {C(n, 6)} p > is, indeed, a polynomial function in the variable 9. Using the fact 
V e S(R d ) and the bounds in Step 3, we have 



3-1 3-1 
i=l i=l 



dOidXidXj < +oo. 



Therefore, by 8.7, Step 4, and dominated convergence theorem, we arrive at the desired result 
and this completes the proof of part (b). 



9. Explicit form of some coefficients. 

In this section, we compute some coefficients explicitly and again show their finiteness by ap- 
plying some basic inequalities arising from Lemmas 6.1 and 7.1. 

For n = and any j > 1, we obtain by applying iterated times (3.1) that 



Cd, a E 


'g-d/2 

1<~2 


(2tt)< 


j+l)d 



1 Ml 



£ °i) ft V(9 l )d9 t d\ l d\ J 



i=l i=l 



v : > +1 (9)de, 



and Cd a — -rJr — • 



(i + i)! 

where we have also used that (47r) d / 2 p^ ) (0) = E \ S l d J 2 

The following Lemma will be useful to prove that the constants appearing in Corollary (1.1) are 
strictly positive. Part of the following proof can be found in [15]. 

Lemma 9.1. Given < a < 2, there exists N a > 1 such that 

l ~l ^ < P(l < S lt9 < N a ), 



where v a = (2 — a)a 2 -"2 2 -° . 

Proof. Let a > be fixed and observe that Si t <* < a if and only if e~ ASl, 7 > e~ Aa , for any A > 0. 
Therefore, Chebyshev inequality tells us that 



P(Si a < a) < inf e 



aA-A 2 



On the other hand, lim < n) = 1 implies that given e > 0, there exists a positive 

n— -f+oo '2 

integer N such that 

1 - P(Si,!>l < n) < e, for all n > N. 



It follows then that for e 



l-e~ v ° 



there exists N a > 1 such that 



1 - P(5i )f < N a ) < 1 | - or P(5i - < ^Va) > 1+ ^ - 
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Now use above facts with a = 1, to obtain that 

P(KS lt9 <N a ) = P{S ll9 < N a )-P{S lt9 < 1) 

1 + e-' u = 



> 



2 

1-e" 



as desired. 



□ 



Remark 9.1. Before proceeding, we give an explicit expression for N a when a = 1. Observe 
that v± = —4. TTie 1/2-subordinator S can be expressed as the first hitting time for the standard 
one- dimensional Brownian motion {Wt} t>0 . More precisely, 



S t = inf js > : W s = -j= 
It is also known (See [1, pp 23-24] for details) that its density is given by 

,(1/2) 



^ (s) = _J_ s -3/2 e - t 74 S> 



2^ 

Therefore, it is not difficult to see that for any N > 1. 



P(l < Si i < N) = 



We can take the 



1 



-3/2 e -l/ 4sds > 



-1/4 /-JV 



2^ J 1 



= -1/4 



3 / 2 ds = L^(i_tv-i/2). 



e- 1 / 4 i/9 1-e" 1 / 4 



equivalently, Nx = jl — -^(e 1 / 4 — 1) j which is approximately 1.786. 



Let us now consider the case n = 1 and j = 2 . We recall that 
(9-1) ^_ (Ai _ A2) + 5^ 1 _ A2 = 5 1 , 

provided < A2 < Ai < 1. In addition, S r *_^ i _^ a \ and S" Al _ A2 are independent random variables. 
Then, it follows by Lemma 6.1 that 



(9.2) 



C {a) {V) = C d,* K i( d ' a ) AV,V) = Cd ^ Kl ^ a ) 



{2n) d x ' ' (2ir) d 
where we have replaced Si by the left hand side of (9.1) to obtain that 



\\7V(6)\ 2 d0, 



KM a) 



1 n \ 



E 



Jo 

1 Mi 



JO 



^l-(Ai-A2)'^A 1 -A2 
( 5 1-(Ai-A 2 ) + ^Ai-Aa)^ 5 

— W W 



d\2d\\ 



(su w + s* w y+i 



dwdXi. 



We now claim that K\(d, a) is both finite and strictly positive when either d — 1 and \ < a < 2, 
or d > 2 and < a < 2 as follows. 

We start with d — 1 and ^ < a < 2. By Lemma 7.1, we obtain in this case that 



< Kx(l,a) < 2- 3/2 E 



c-l/4 
°l,a/2 



1 r X 



Jo 



{(1 — «;)«;} 2a dwd\\. 
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On the other hand, when d > 2 we have 



(9.3) 



Q<Kx{d,a) < 



l 



1 

2 Jo Jo 



E 



\°l-w "r °w) 



dwdXi 



E 



-^l-d/2 
'l,a/2 



where we have used the basic inequality 

ab 



ab 



<ka + bf- d l^ 



(a + 5)^/2+1 (a + by^-^a + b) 2 ~ 2 

valid for all a, b > 0. The expectation in (9.3) is finite for all a since 1 — | < < ^. Therefore 
Kx(d, a) is finite in the cases stated above. 

Next, we prove that K\(d, a) is strictly positive. We note that (2.1) implies that St = t*S\. 
Therefore, we can write 

S{_ w = (1 - w)ix u S* w =wix 2 , 

where X\,X 2 are independent copies of S\. That is, X\ = S\ = X 2 . Thus, for any < w < 1, we 
have 



^1 — w^w 



= E 



{l-w)iw 2 / a X 1 X 2 
{(l~w) 2 / a X 1 +w 2 / a X 2 ) 1 + 

X\X 2 



> (1- w)«w~E 
(1 — w) = w 



> 



(X 1 + X 2 ) i +* 
^^(1 < Si, f < iV«) 2 



T ;KX 1 ,X 2 <iV Q 



{2N a ) 

From Lemma 9.1 and last inequality, we conclude that 

,— v„\2 r 1 r x i 



K\ (d, a) > 



(1-e- 



(1 — id) <*w a dw d\i > 0. 



4(2iV a ) 1+ f Jo Jo 

Similarly, for either all a G (1,2) and d > 2, or for d > 4 and a € (0,2), we have 

\AV(6)\ 2 dd, 



pWrtrt _ Cd, a K 2 (d,a 
°2,2 l K i — 



where 



# 2 (d,a) 



1 /■ A 



2(2?r) d 

/ Q* \2 

Wl-(Ai-A 2 )°Ai-A 2 / 

2+i 



^ L( 5 1-(A!-A 2 ) +^A 1 -A 2 ) 

By applying the same argument as above, we have the following 
(i) a E (1, 2) and d > 3 or for d > 4 and a £ (0, 2), we obtain 



d\ 2 d\i. 



-v a \2 r 1 r x 



< 



(l-e~"°) 
4(27V a ) 2 +§ 7 Jo 



(1 — w) w dw d\\ < K 2 (d 1 a) < 



S 



2-d/2 



12 



(ii) a € (1, 2) and d = 2, we have the same lower bound as in (i), but by Lemma 7.2 we obtain 



fsT 2 (2,a) < 2"' j £; 



5" 



1/2 

l,ct/2 



JO 



{«;(! — w)} a dwdXi 
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Likewise, it is not hard to prove that 



where 

pi p\\ p\2 

K 3 (d,a) = E 



JO JO 



J l-(Ai-A 3 ) °Ai-A 2 + l-(Ai-A 3 ) A2-A3 + ^At-Aa J A 2 -A 3 

rs* rx n + s: x +sx x ) 1+d/2 

\ 1— (Ai — A3) Al — A2 A2 — A3 / 



d\3d\2d\i 



is positive and finite provided that either d > 2 and 0<a<2or<i=l and i < a < 2. 

For the rest of the paper, we will use the following notation for the constants given above, 

C d . a K 3 (d,a) ., C d , a K 2 {d,a) 

Md >« = and Nd - a = w • 

Remark 9.2. Based on the computations in [6] and part (b) of Theorem 1.1, we have under the 
conditions stated above that 

limAfd.a = 7^7: and limMd - = 777- 

OC\2 120 0£\2 12 

10. Proof of Corollary 1.1. 
The proof uses a combination of Theorem 1.1 and Lemma 7.2. 

Case M=l. 

1) When d > 2, we invoke Theorem 1.1 with J — 3. We have in this case 1 — | < < 
Therefore, we can consider any a on the right hand side of next expression, 



(10.1) 



Tr(e tH - e + ( f v{g)dg _t* f y2( ^ + g /" y 3 (0)d0 = ^i Q) (i)#(«) 



^ a) (o) jR d 2 At'* 3! - 

where ^4 (1) = min {4, 2 + — }. Hence, (iii) follows by noticing that 



2 2 
2+-<4<3+-, when a e (1,2) 
a a 

2 2 
4<2 + -<3+-, when a € (0, 11. 
a a 

2) For the case d = 1, we use Lemma (7.2) which guarantees that (10.1) is still true for 
i < a < 1. Thus (i) holds. 

Case M=2. 

We apply Theorem (1.1) with J = 4 and 2 - | < f to obtain 

rr(e ~ t t ) rr tH " ) +* / vy)d0-~[ v\e)dd + c$(v)t 2 ^ 
(10.2) +^ / v 3 (e)de-c i 1 a i(v)t 3+ ^ ~^ f v 4 (e)de + c[ a }(v)t 4+ i = t^^R^(t), 

where ^ q) (2) = min {5, 2 + ^}. 
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1) (iv), (v) and (vi) when d > 4. 

We have 2 — | < < % and then any a can be considered. For a £ (0, 1], we have 
2 < — , which implies that all the powers of t containing a are larger than 5 except 2 + — . 
Comparing 2 + — with 5 yields (v) and (vi). Notice that (iv) also follows since all the 
power of t containing a in (10.2) are larger than 4 (simply use the fact that 1 < — ) except 
2 + — which is less than 4 whenever 1 < a < 2. 

a 

In addition, by Lemma 7.2 we obtain that (10.2) remains true in the following cases 

2) When d = 3 and \ < a < 2. Hence, (iv) and (v) holds by part 1). 

3) When d=2, < a < 1. Then, (iv) also holds for d = 2. 

4) When d = 1, § < a < 2. Thus, (ii) holds. 

This covers all the cases and completes the proof of Corollary 1.1. 



11. Explicit expansion for a = 2/k, k > 2 integer and a close to 2 

In this section we want to provide to the reader a better insight of our main theorem by finding 
a expansion formula of the trace when a — -| with k > 2 an integer and for values of a near 2, of 
course, under the condition that 2M — d < a, which is equivalent to 2M — d < when < a < 1. 
We also give examples as an application to the results given below. We refer to the reader to the 
end of §9 for the definition of the constants Cd. a , M-d.a and Afd,a- 

Theorem 11.1. Let a = 1 and 2M - d < 0. Then, for any 2 < J < 2M + 1, 

J ( \ 



(11.1) 



Tr(e 



t / v(e)de-^2t l Yl (~ l ) n+3C nW) 



1=2 2n+j=l, 
\ J>2 



0(i J+1 ), 



as 1 1 0. 



Proof. We apply Theorem 1.2 with a = 1 and J+l < 2M+2, so that ^Jyj^M) =min{J+l,2Af + 2} 
J + 1. Therefore, we obtain 



Tr{e- H ^ - e~ Hlt ) 



+ t 



V(d)dd 



E 

j+2n<J, 
2<j<J, 0<n<M-l 



(-l) n+ ic£](V)t 2n+j = 0(t J+1 ). 



as t I 0. 

The following argument shows that under the conditions given above we are not excluding any 
n and j such that 2n + j < J. We observe that 2n + j = I for some I E {2, J} if and only if 



1-3 



Then the larger n can be is 



.7-2 



But J < (2M + 1), which yields ^ < M 



h. Since 



n is a positive integer, we conclude the larger that n can be is in fact M — 1. This implies that 



E 

3+2n< J, 

2<j<J, 0<n<M-l 



1=2 



. 2n+j=i 

\ i>2 



/ 



□ 
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Example 11.1. When M = 2 and d > 4, we have 2-2-d<0<l . Then, Theorem 11.1 holds 
for any 2 < J < 5. Therefore, for the particular case J — 5, we obtain 



Tr(e- H ^ - e- Hlt ) 



+1 j V{9)d9- t - I V 2 (9),W 



4L v3mM -UL r 



(0)d0 + 4\C dA / \WV(9)\ 2 d9 



(11.2) 



5! 



V 5 (9)d9 - 5!M d ,i / V{9) \VV{9)\ 2 d9 = 0{t 6 ), 



as t 4- 0. Notice that part of this expansion is obtained by applying v) of Corollary 1.1 to the specific 
case a = 1. 

By mimicking the proof for the case a = 1, we conclude that 

Theorem 11.2. Lei ct = \ with k > 3 a positive integer. Assume also 2 AT — d < 0. Then, for 
any 2 < J < 1 + Mfc we /iai>e 



(11.3) 
as fj-0- 



Trie 



-Hvt 



-H 2 t. 



( 



P <*>(o) 



t / v(0)d0-]T)t' 



\ 



2 (-lr+^'ao 



kn-\-j—L, 
\ J">2 



0(4 



Example 11.2. Consider k — 3, M = 2 and d > 4. Then, our main theorem holds for 2 < J < 7. 
XTie particular case J — 5 yields 

-h 2 t s 



Tr(e 



-H v t 



3 ) 



P^(0) 



4 d 



t / y(0)d0 _ _ / v^(0)d0 + ^ / ^(0)d0 



3! 



(11.4) 
as 1 10. 



4! 



y 4 (6»)d6» 



5! 



7 b (f)df + 5!£ d; 2 / |VV(fl)rtW = 0{t b ), 



Let us consider for J > 2 the following J — 1 x J — 1 matrix which contains all the power of t with 
the form — + j that may appear in the expansion of the trace, Aj(a) = (a r , s ) with 1 < r < J — 1 
and 

r-s + 2+ |(s-l) if r < s, 
otherwise. 
In this matrix, n = r — s + 2 and j = s — 1. Observe that n + j = r + 1. Thus, 

Example 11.3. For a = l, 

\ ( 2 ■ • • • 
'34.. . 

4 5 6. . 

5 6 7 8. 

6 7 8 9 10 



A 6 (a) 



/ 2 












3 


2^ 


- 1 


2_ 

a 






4 


3^ 


- 1 


2 

a 


2- 


-2 


5 


4^ 


- 1 


2_ 

a 


3- 


-2 


6 


5H 


- 1 


2 

a 


4- 


-2 



^ 2 + 3- ^ 

a a 
2 o 1 o 2 



V 
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We have set the two matrices together to match entry by entry. As an example, we conclude that 
there are two coefficients related to t 5 . Namely, CqI(V) and C^l(V). The reader can verify this 
conclusion from 11.2. 



Likewise, for a = I we obtain 



( 2 










3 


5 








4 


G 


8 






5 


7 


9 


11 




6 


8 


10 


12 


14 . 


V 7 


9 


11 


13 


15 17 J 



A 7 (2/3) = 



We can deduce then that the next two terms in the expansion given in 11. 4 are 

{{-i)*c$(v) + (-i) 3+1 cgV)} - 1 7 {(-i) 7 c$(v) + (-i) i+1 c£>(v)} . 

We point out that for any < a < 2, we always have 2 < 3 < 2 + — which implies that the 
influence of the a in the expansion of the trace is expected to be seen in some place after the term 

Notice that for every J > 2 we have 



(11.5) 



A/(2) 



( 2 






\ 


3 


3 






4 


4 


4 . 




5 


5 


5 5 




V J 


J 


J J . 





which says, for example, that in the expansion (1.2) there are three coefficients associated with 



t 4 . Namely, C$(V), C$(V) and CfflV). 



(2), 



Theorem 11.3. Assume J > 4 and 2(J — 2) — d < 1. Then, for all a € ( 2( /_ 3) ,2) we have 



Tr(e 



-tH v 



tf V(9)d9-J2 E (-l) n+J C { n a ](V)t^ =(D(t J ), 

JR d , o _ , , , 



1=2 n+j=l 
3>2 



as t J, 0. 

fVoo/. In Theorem 1.1 we take J = M + 2 > 4 so that 



*£*(M) = $g x (J - 2) - min ( J + 1, 2 + 2(J 2) f > J 



since — > 1. Now we want to choose a such that a TtT < a r +i,i for r G {2, J — 2}. This last 

,7-2 



condition is equivalent to < a and implies that all the entries of Aj(a) are increasing for 
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these ct's. In other words, we have the following arrangement 

2 

(11.6) <3<2 + - 

a 

2 2 
<4<3 + -<2 + 2- - 
a a 



< J - 1 < (J - 2) + -<...< 2 + (J - 3) • - < J. 

a a 



Then by Theorem 1.1, 

Tr(e- tHv - e~ tH ") 



+ t f V(8)d8- Yl (-^) n+3 C^(V)t 



2<j"<J-l, 0<n<J-3 



(11.7) (-t)*^ 1 ( J - 2 ) J R. 7+1 (t)+ £ 



2<j<J-l, 0<n<J-3 



We observe that the right hand side of (11.7) is 0(t J ) as t I 0, due to $,/+i(J — 2) > J. On 
the other hand, it is easy to see by the definition of Aj(a) that the only powers of t satisfying 
— + j < J are those in the arrangement given in (11.6). Therefore, due to this arrangement we 
can also rewrite the third term in the left hand side of (11.7) as follows 



j-i 



(11.8) ]T (-l)»+icg(V0^+i = £ £ (-ir^C^(V)t^. 

2<j<J-l, 0<n<J-3 i^ 2 



□ 



Example 11.4. We take J — A and d > 3 so that 2 • 2 — d < 1. Then, according to last theorem 
for all a G (1, 2) we have 

>-+t / V(0)dO- - / V 2 (9)d9 



t 3 

+ "3! 



("Vrtt .Ld 2! 



\[ V 3 (8)d8 + £ d , a t 2 +^ [ \VV(8)\ 2 =0(t 4 ), 



as t J, 0. Notice that this result is already given in Corollary 1.1. 

Let us now consider J = 5 and d > 5. Then, for all a € (§i 2) we obtain 

Tr ( e -Wv _ e -tH a) , f2 , / t 3 



■ * / V(8)d9- t -[ V 2 (9)d9+(^f V 3 (8)d9 + C d , a t 2 +^ [ \VV(9)\ 2 ) 
V 4 (8)d8 + M dM t 3+ i [ V(9)\VV(9)\ 2 d9 + M d . a t 2+ ^ [ \AV(9)\ 2 dff) = 0(t 5 ), 



*T(o) 

4! 



as i | 0. 
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12. Extension to a-RELATiviSTic processes. 



In this section we describe how to use the tools developed in the preceding sections to obtain 
an asymptotic expansion for relativistic stable processes. 

Let < a < 2 and m > 0, we consider the function (f) m (\) = {A + to 2 /"}"^ 2 — to. As the 
case Tti — was already studied above, we assume for the rest of the paper that m > 0. The 
function <f> m is known to be a Bernstein-function, that is (— l) n (j)m\^) < for every n £ N and 
A > 0. Therefore, there exists a unique subordinator {S t , m } t>0 (see [ : , pp 89]) such that its 
Laplace transform is given by 

E[e~ XSt -™] = e -*^( A ). 
It is easy to see that the transition density rfP(s) of St, m satisfies the scaling 

e mt r,l a/2) ( S )e- m2/as , 



TO 

which implies that for every — oo < rj < a/2, 
(12.1) 0<E[Sl] m ] = e mt E 



< e mt t^E 



J 1,q/2 



< OO, 



where S^st is an a/2-subordinator. 

The a-stable relativistic process in M. d is defined as the subordinated Brownian motion process 
Z™ = i?2S t>m , whose characteristic function is given by 

(12.2) E\e- l< ^>] =e -*^(l«l 2 ). 

If p£ Q ' m ' (a;) denote the transition density of Z™ then exactly as in the preceding sections we 
have 

pt m \x)=E[p$Jx) . 
Consequently, because of (12.1), it follows that 

(12.3) t d ' a p ( t a ' m \0) = (ATt)- d ' 2 E[s- d t lf 

= {Air)- d/2 e mt E 

which implies by monotone convergence theorem that 
and for all < t < 1, 



5 -d/2 e -(tm) 2/Q S 1 _ 



t|0 



(12.4) 



t d / a p[ a ' m) (0) > {^)- d ' 2 e-' m E 



si 



-d/2 



We point out that in order to obtain the constants in the asymptotic expansion in Theorem 
(1.1), the self-similarity of the subordinator St, a, namely that t 2 l a S c ^ = ^ct.f and the scaling 
property of p[ a \x), were strongly used. But, these two properties are not satisfied by St, m an d 

(a,m) 



p\ ' {x). Indeed, we have 
t 2/a S CA 



v 



S' 



and p ( ^ m \x)^m d / a pt t 1) (m 1 / a x). 



Despite of this, we can still obtain a version of Theorem 1.1 where the coefficients C„j (V) are 
replaced by time dependent functions. 
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We now define the following linear operators H a ^ m = cf> m (—A) and Hy — H a . m + V. In [16] it 
is proved that for any Bernstein function the Feymann Kac formula holds. As a result, the heat 
kernel of Hy can also be written as 



(a,m) / 



<J x,y 



-ftV(Z™)d 8 



We now proceed to mention which of the results in the above sections still hold true for these 
operators. It is easy to check that all the results in section §3 and §4 are still true when we replace 
ICT by m (IC| 2 ), Pt by p { t a ' m) and using the fact that m (-A)/(f) = <£ m (|£| 2 )/(£). Furthermore, 
by replacing the random variables S\ by S\^m and by tm , we see that all results in §5 also 
remain true with the following modifications. Set, 



^ ( a,(,^r (l " i)Wltf, ^ ,At - Vi 



and 



3-1 



Lj(tm,X,9) = J2 S L-\ k +i,t m \lk\ 2 



k=l 



Si, 



3-1 
k=l 



where A = (Ai, Xj) satisfies (5.5) and 



n,— 



Then 



We recall that we cannot express the right hand side of the last equality in terms of p\ a (0) 
because this density fails to satisfy the scaling property. We now state the analogue of Proposition 
(3.1) for the new linear operators defined above. 

Proposition 12.1. Given J > 2, we have 



Tre 



-tH v 



-,-tH a 



Pi 



(a,m) 



(0) 



= -t 



V{0)d0 



/•Xi 










f E 


10 


'Jo J 





o-d/2 -t > / a L i (tm,\,0) 
°l,tm e 



+ § (27T)W«j>J a,mj (0) Jo Jo 

3-1 3-1 

x V(-J20i) YlVi^deidXidXj + (-t) J+1 R J+1 (t), 

i=l i=l 

where Rj+i(t) is a bounded function on (0, 1). 

From this lemma, it is clear that with the Taylor-expansion of the exponential function we 
obtain time depending functions instead of coefficients as in §6. Now, it is easy to see by using the 
expansion (6.1) of order 2 that when 1 < the following equality holds 



E 



si 



■d/2 -t 2/a Lj{tm,\,8) 



t*/«pj a ' w °(0) 



(Air) d/2 -t 2 / a R ld (t), 



where 



E 



R 1J (t) 



S-£L j (tm,X,e)e-^ 
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for some nonnegative random function f3\ = (3i(tm, X,9) and this satisfies, according to (12.4) and 
(12.1), the estimate 

( 4n y/2 e 2m E \ S 7 d J 2+1 ] 3-1 
& [ b l,m \ fe=l 

for all < t < 1. With this last estimate, by mimicking the proof of Theorem (1.1) ,we have that 



(12.5) 



Tr{e- Hvt - e- Hmt ) 



Pi 



(a,m) 



(0) 



+ t / V{9)d9~ 



2! 



V 2 (9)d9 



3! 



V\6)d6 = t** ) VR { r ) {t), 



where 3>£*^(1) = min {4, 2 + and R^it) is a bounded function of t over (0, 1). We can also 
redo the proof of Lemma 7.2 to make sure that (12.5) remains true when d = 1, M = 1, and 
i < a < 2. For this, we only require the following inequality when < t < 1, that follows again 
from (12.1), 



E 



S, 



M/2-d/i 
Ltrn 



[ M/2-d/4 E 



S 



M/2-d/i 



l,ltm 



< l M/2-d/i e lm E 



M/2-d/i 



From this we conclude 
Corollary 12.1. 

Tr(e~ Hvt - 



Pi 



(a,' 



'(0) 



+ t I V(6)d6-y I V 2 (9)d9 + - j V 3 (9)d9 



3! 



0{t 2 +«), if a € (1,2) and d>l, 
0{t 4 ), ifae (§, 1] and d>l, 
C(i 4 ), tfae (0, i] andd>2, 

as 1 1 . 

This extends the result in Bahuelos and Yildirim [3] where the second term is computed. 

13. Extension to mixed-stable processes. 

In this section we assume that 0<a</3<2,a>0 and only extend Theorem 1.1 for the case 
1 < ^±zL, In addition, we analyse the cases d — 1, 1 < a < 2 and d = 2, < a < (3 < 2. We 
point out that the constants C a ,d,a,/3 below may change from line to line. Consider the Bernstein 
function </> a (A) = A 2" + and let {St,a} t>0 be the subordinator associated to 4> a , with Laplace 
transform given by 



E[, 



,-ASt.al _ p -t<j> a (\) 



It is a standard fact that the above subordinator can be written as 



(13.1) 



St,a — Si, 4 



its. 



a- 



where the subordinators in right hand side of the last equality are independent processes. 
Notice that, by independence, we have 

(13.2) St ia =t*Si, 9 + (at)$S lA , 

■2 -L, 2 

and by (13.1), it also follows that 

{l + a?)mm{s t ,f,S t ji} < S t , a < (1 + a?) max {X f , S t | } , 
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which implies for any rj < #, 



(13.3) 



0<E[Sl a ] <2(l + a?)"max{.E 
<2(1 + a?) 7 ' max jt^? 



SI, 



,E 



E 



2>7 

ATE 



1 2- 



< 00. 



The (i-dimensional process Y" t a = i?2S t a is known as the mixed (a, /3)-stable process, with 
characteristic function given by 



E e -i<t,v t a > 



-t<MI?l 2 



Properties for the heat kernels and Green functions for these processes have been studied by 
several people, see for example Chen and Kumagai [ I], Chen, Kim and Song [9] or Jakubowski 
and Szczypkowski [17]. Using some of these results, Banuelos and Yildirim [3] were able to compute 
the second coefficient. In this action we show how to go further than the results in [3] using the 
techniques of this paper. 

Denote the transition density of Y t a by p^ (x) , then as before we have 



p[ a \x)=E\pf ta (x) 



Therefore, by (13.2), the following holds 



Pi 



(a) 



(0) = {4ir)- d / 2 E 
= {Att)-^ 2 E 



n-d/2 



t*Si, 9 + {at)fS lA 



-d/2 



It is also proved in [9], that there exists a constant Cd, a ,f3 such that 



where 



/f(x) = min J min (ai)"^} , — ^ + 



at 



In particular, it follows that there exist a constant C a ,d,a,/3 > such that 



(13.4) P?\Q)>C aM r d ! fs , 
provided < t < min jl, a' J -= j. 
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Set 



F {a >(tX,^6) = e fc=i '=1 



A (i, A) = ^5r_ (Al _ Aj)>f + (ai)^*_ (Ai _^ |: 



3-1 



4 Q) (t,A,0)=^A fc (i,A)| 7fc f 



3-1 

E A fc (t,A) 7fc 
fe=l 



fe=l 

Then, we obtain by using the tools of §5 that, 
(13.5) f F ia \t\,^9)e~ tx ^^ ) d£,='K- d/2 E 
We also observe that 



3-1 



(13.6) 



With this equality, we can imitate the proof of Lemma (6.1) along with the Taylor expansion to 
arrive at 



E 



(13.7) 
with 



< 



-d/2+1' 



3-1 



^ a) (o) 



fc=i 



3-1 



(13.8) 



fc=i 



where we have used (13.4) and (13.3). 

Let us set H a = (f) a (—A) and Hy = H a + V. With this notation and the above bounds we 
obtain (by mimicking the proof of Theorem (1.1)) the following result. 

Theorem 13.1. Assume d=l and 1 < a < (3 < 2. Then, 



Tr(e 



— Hyt ~ — tH, 



' +t I V(9)d9- t - I V 2 (9)de+ t - I V 3 (6)d9 = o{t 2+ f). 



as 1 1 . 
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Theorem 13.2. Assume d>2,0<a</3<2 and 2 + § - d ^~ - > 0. Then, 



p^iO) Jm" 2! 7 Rd 3! 

/or some bounded function R{t), t G ft), min |a^° , l|j . Here, 



$ d (a, /3) = min <j 4, 2 + - - d [ - - \ 
a \a p 



As an example, observe that 

$2 (<*,#) = min ^4, 2 + 

This gives 

Corollary 13.1. For d=2, we have 

T ( -H v t _ -tH a \ r t 2 p t 2 

[ (ah ' l + t v(d)dd- t - v\e)de + t -f v\6)de 

\ a > (\\ 7b2 2! ./„2 3! 



C(i 2+ 7), if 1 < f3 and a < (3, 
0{t 4 ), if < a< & < 1 . 



Finally, we remark that the condition 2 + ^ > d — comes from (13.8) in order to obtain 
the boundedness of R(t) over the interval stated above. 
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